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ABSTRACT 

Taking advantage of the very precise de Jager et al. optical white dwarf orbit and spin 
ephemerides; ASCA, XMM-Newton, and Chandra X-ray observations spread over 10 
yrs; and a cumulative 27 yr baseline, we have found that in recent years the white 
dwarf in AE Aqr is spinning down at a rate that is slightly faster than predicted by 
the de Jager et al. spin ephemeris. At the present time, the observed period evolution 
is consistent with either a cubic term in the spin ephemeris with P = 3.46(56) x 
10^^^ d~^, which is inconsistent in sign and magnitude with magnetic-dipole radiation 
losses, or an additional quadratic term with P — 2.0(1.0) x 10^^'^ d d^^, which is 
consistent with a modest increase in the accretion torques spinning down the white 
dwarf. Regular monitoring, in the optical, ultraviolet, and/or X-rays, is required to 
track the evolution of the spin period of the white dwarf in AE Aqr. 

Key words: binaries: close - stars: individual (AE Aqr) - novae, cataclysmic variables 
- stars: rotation ~ white dwarfs. 



1 INTRODUCTION 



AE Aqr is a bright {V ~ 11), nova-like cataclysmic bi- 
nary consisting of a magnetic white dwarf primary and a 
K4-5 V secondary with a long 9.88 hr orbital period and 
the shortest kno wn white dwarf spin period P — 33.08 s 
JPattersonl Il979r) . Although originally classifi ed and inter - 
preted as a disk-accreting DQ Her star (.PattersonI Il994) . 
AE Aqr displays a number of unusual features that are not 
naturally explained by this model. First, violent flaring ac- 
tivity is observed in the radio, optical, ultraviolet. X-ray, 
and TeV 7-rays. Second, the Balmer emission lines are sin- 
gle peaked and produce Doppler tomograms that are not 
consistent with those of an accretion disk. Third, the white 
dwarf is spinning down at a rate P = 5.64 x 10~^* s s^^ 
(|de_Jager et al. 1994, hereafter, de Jager). Although this 
corresponds to the small rate of change of 1.78 ns yr~^, 
AE Aqr's spin-down is typically characterized as "rapid" 
because the characteristic time P/P ~ 2 x 10^ yr, which is 
short compared to the lifetime of the binary, and because 
the spin-down luminosity Lsd = —IftO, ~ 1 x 10''^ erg s~^ 
(where I « 0.2Mwdi?wd ~ 2 x 10^° g cm^ is the white dwarf 
moment of inertia, M„d and -Rwd are the white dwarf mass 
and radius, respectively, Q = 27r/P, and tl = —2-kP/P^), 
which exceeds the secondary's thermonuclear luminosity by 
an order of magnitude and the accretion luminosity by two 
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orders of magnitude. Given this, AE Aqr could be thought of 
as being powered primarily by the ultimate in clean energy 
sources: a flywheel. 

Because of its unique properties and variable emis- 
sion across the electromagnetic spectrum, AE Aqr has 
been the subject of numerous studies, including an inten- 
sive multiwavelength observing campaign in 1993 October 
JCasares et al.l ll996. and the series of papers in ASPC 85). 
Based on these studies, AE Aqr is now widely believed to be 
a former supersoft X-ray bin ary Ij Schenker et al.ll 2002ll an d 
current magnetic propeller iWvnn. King, fc Hornd 119971) . 
with most of the mass lost by the secondary being flung 
out of the binary by the magnetic field of the rapidly ro- 
tating white dwarf. These models explain many of AE Aqr's 
unique characteristics, including the fast spin rate and rapid 
secular spin-down rate of the white dwarf, the anomalous 
spectr al type of the secondary, the ano malous C to N abun- 
dance JMauche. Lee, fc Kallm an 1997' ), the absence of sig - 
natures of an accretion disk ||Welsl^Jlom«^^^om«i^^9^), 
the violent flaring activity JPearson. Home, fc Skidmord 
|2003^, and the origin of the radio and TeV 7-ray emission 
(KuiiDcrs et al. 1997; Mcintics & Venter 2003). 

To build on this observational and theoretical work, 
while taking advantage of a number of improvements in ob- 
serving capabilities, during 2005 August 28-September 2 a 
group of professional and amateur astronomers conducted 
a campaign of multiwavelength (radio, optical, ultraviolet. 
X-ray, and TeV 7-ray) observations of AE Aqr. Analyses of 
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these data are ongoing, but here we present a fundamental 
result - the spin period of the white dwarf - that relies solely 
on photometric data from the Chandra X-ray Observatory 
and archival data from ASCA and XMM -Newton. 



2 OBSERVATIONS AND ANALYSIS 

AE Aqr was observed by ASCA with the Gas Imaging Spec- 
trometer (GIS) detectors and the Solid-State Imaging Spec- 
trometer (SIS) detectors beginning on 1995 October 14 at 
00:16 UT for 82 ks (Sequence Number 33005000), by XMM 
with the European Photon Imaging Collaboration (EPIC) 
pn CCD detector beginning on 2001 November 7 at 23:47 
UT for 14 ks (ObsID 0111180201), and by Chandra with the 
High-Energy Transmission Grating (HETG) and the Ad- 
vanced CCD Imaging Spectrometer (ACIS) detector begin- 
ning on 2005 August 30 at 06:37 UT for 81 ks (ObsID 5431). 
Th e ASCA a n d XMM data have been previous l y disc ussed 
bv lEracleou^ ^1999^ : IChoi. Dotani. fc Agrawall il999D : Os- 
borne [2002, unpublished presentation at the Third Mag- 
netic Cata clysmic Variable Workshop (lAU Coll. 190)]; and 
lltohet al.U2006l) . 

The ASCA and XMM data and the Chandra data 
were extracted from the High Energy Astrophysics Science 
Archive Research Center and the Chandra Data Archive, 
respectively. The data files used for subsequent analysis 
were the ASCA GIS and SIS screened bright mode event 
files, the XMM EPIC pn pipeline processed event file, and 
the Chandra level 2 pipeline processed event file. The data 
in these files were manipulated as follows. First, all times 
were converted from spacecraft time to Terrestrial Time 
(TT) and corrected to the solar system barycenter using 
the ASCA FTOOLS^ v6.0 tool timeconv, the XMM Sci- 
ence Analysis Software (SAS'^) v6.5.0 tool barycen, and the 
Chandra Interactive Analysis of Observations (CIAO^) v3.2 
tool axbary. The barycentric corrections produced by these 
different software packages were checked against those of 
the Interactive Data Language (IDL) procedure barycen* 
(which does not by default account for the varying light 
travel time between the satellite and the center of Earth), 
and were found to agree within 40, 210, and 200 ms for 
ASCA, XMM, and Chandra, respectively, consistent with 
the size of each satellite's orbit. Second, source and back- 
ground events were extracted from the event files using 
custom IDL software. For Chandra, events were collected 
from the zeroth order image and the ± first order dispersed 
spectrum using the region masks in the level 2 pha file. 
Third, event times t were converted to white dwarf orbit 
and spin phases via the relations (j>orb ~ floih{t — To) and 
<^spin = 0,o{t - Tmax) + ^^^(i ~ Tmax)^, respectively, where 
r^orb = 27r/Porb, ^Q = 2tv/P33, aud n = -27rP33/P33 and 
To, Porb, rmax(BJD) ^ 2445172, P33 = 0.00038283263840 d, 
and P33 = 5.642 x lO^^** d d~^ are the white dwarf orbit and 
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Figure 1. (a) Chandra spin phase folded count rate light curve 
of AE Aqr {filled circles with error bars), best-fitting cosine [solid 
curve), and mean count rate A [dotted line), [b) Spin phase offset 
as a function of orbit phase [filled circles with error bars), best- 
fitting cosine [solid curve), mean phase offset (dotted line), and 
the predicted spin phase offset variation for the de Jager white 
dwarf orbit ephemeris and a pulse time delay of 2 s [dotted curve). 



spin ephemeris constants from Table 4 of de Jager. ^ Fourth, 
filters were applied to restrict attention to events from two 
orbital cycles for ASCA and Chandra and | of an orbital 
cycle centered on 0orb = 0.25 for XMM. The resulting range 
of orbit phases and observation dates are listed in Table 1. 
Fifth, background-subtracted spin phase folded count rate 
light curves were calculated and fit with the cosine function 



is the pulse semiamplitude, and 4>o is the phase offset (ex- 
pressed in cycles (f>o/2n in what follows). Given that the op- 
tical, ultraviolet, and X-ray spin pulses are aligned in phase 
(Patter son et al. 1980 ; Eracleous ct al. 1995) , (j)o should be 
equal to zero to within the uncertainty of the de Jager spin 
ephemeris, if the ephemeris remains valid at the times of the 
ASCA, XMM, and Chandra observations. 



1 Available at|http://heasarc. gsfc.nasa. gov/docs/software/lheasoft/ftoots/ 1ote that the expression in §6 of de Jager for the times of spin 

2 Available at http://xmm.vilspa.esa.es 7s as/| maxima formally introduces a term proportional to [t — Tmax)^ 

^ Available a t ,http: //cxc. harvard. edu/ciao/| in the expression for (/ispin, but it amounts to an insignificant 

* Available at|http://astro. uni-tuebingen.de/software/idl/aitlib/astro/ [2 x 10~^ cycles at the time of the Chandra observation. 
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Table 1. Log of observations. 



Satellite 



Orbit Phase 

<t>ovh (cycles) 
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Date 
[BJD{TT) - 2400000] 



ASCA 11738.8000-11740.8000 

XMM 17124.0625-17124.4375 

Chandra 20504.0000-20506.0000 



50004.62127-50005.44459 
52221.49479-52221.64916 
53612.86503-53613.68834 



2.1 Chandra results, ignoring the pulse time 
delays 

We begin the analysis by considering the Chandra data, 
which covers two full binary orbits of AE Aqr without in- 
terruption by Earth occultations or detector shutdowns. As 
shown in the top panel of Figure 1, the Chandra background- 
subtracted spin phase folded count rate light curve of AE 
Aqr is reasonably weU fit {xl = 19.1/17 = 1.12) by the 
cosine function with a mean count rate A = 0.313 ± 
0.002 counts s~^, a relative amplitude B/A = 15 ± 1 per 
cent, and a phase offset (f>o = 0.232 ± 0.011 cycles (through- 
out the paper, errors are la or 68 per cent confidence for 
1 degree of freedom); the phase offset differs from the de 
Jager spin ephemeris by 4.4(t. To check if this result is af- 
fected by intensity variations around the orbit, we repeated 
the above procedure for each of 8 contiguous orbit phase 
intervals centered on (jiorb = [0, 1, 2, . . . , 7]/8. The result- 
ing variation of the spin phase offset with orbit phase is 
shown in the bottom panel of Figure 1 and is well fit (xt = 
4.31/5 = 0.86) by the cosine function A' + B' cos((^orb — i^o), 
with A' = 0.234 ± 0.010 spin cycles, B' = 0.066 ± 0.015 
spin cycles, and (pQ = —0.03 ± 0.03 orbit cycles: the orbit 
phase offset (/jq is consistent with zero, the mean spin phase 
offset A' is consistent with the value derived above for (f>o, 
and the semiamplitude B' corresponds to a pulse time delay 
B' Pzz = 2.17 ± 0.48 s, which is consistent with that mea- 
sured in the o ptical (2.04 ±0.13 s, d e Jager) and ultraviolet 
(1.93 ± 0.03 s. lEracleous et al.lll994h . From this resuh,'' we 
conclude that the source of the pulsating X-rays follows the 
motion of the white dwarf as it orbits around the center of 
mass in AE Aqr. 



2.2 ASCA, XMM, and Chandra results, 
accounting for the pulse time delays 

Given these results, we analyzed the ASCA, XMM, and 
Chandra data for AE Aqr assuming that in each case the 
X-ray source follows the motion of the white dwarf. Specifi- 
cally, we corrected the barycentric event times by 2 cos{4>ovh) 
s before calculating the spin phases, the spin phase folded 
fight curves, and the cosine fit parameters. The pulse time 
delay correction allows us to use all the data for each obser- 
vation, without regard to gaps due to Earth occultations or 



S iReinsch et alJ ^1995^ previously reported that the ROSAT X- 
ray pulse time delays vary with the orbital period and are in phase 
with the optical/ultraviolet pulse time delays, but they did not 
provide details. 



detector shutdowns {ASCA) or to the varying source inten- 
sity (all three observations). With the pulse time delay cor- 
rections and the previous orbit phase filters, the spin phase 
folded count rate light curves and the best-fitting cosines are 
as shown in Figure 2 and the best-fitting cosine parameters 
and the Barycentric Julian Dates of the X-ray pulse maxima 
are as listed in Table 2. Figure 3 plots the phase offsets ver- 
sus time, demonstrating that the observed phases O are di- 
verging from the calculated phases C assuming the de Jager 
spin ephemeris: accounting for the error in the ephemeris, 
the two differ by 1.7cr, 3.9cr, and 4.4cr for ASCA, XMM, and 
Chandra, respectively. Apparently, the white dwarf in AE 
Aqr is slowly but progressively "coming off the rails" of the 
de Jager spin ephemeris. 



3 UPDATING THE AE AQR WHITE DWARF 
SPIN EPHEMERIS 

The trend displayed in Figure 3 of the difference between the 
observed and calculated AE Aqr white dwarf spin phases is 
steeper than can be explained by an increase in either the 
value of the P33 parameter or the size of the parameter un- 
certainties in the de Jager spin ephemeris. Given the limited 
baseline and the limited amount of X-ray data, the observed 
trend is consistent with the addition to the ephemeris of 
either a P33 term centered on Tmax or an additional P33 
term centered on a late date Tniax- In the first case, the 
Taylor expansion of the spin frequency Q.{t) = f2o + f2(i — 
Tmax) + \^{t — Tmax)^ results in a cubic term in the expan- 

Sion for 0spm = no(t-rmax) + ^n(t-r.^ax)^ + ^n(t-rn,ax)^ 

jManchester fc Taylor! [l977ll . so O - C = -\^{t - T^^^f 
and the optical and X-ray data are consistent with Cl — 
-1.48(24) X 10"" d"^ or P33 = 3.46(56) x 10"^^ d"\ In 
the second case, O — C = —^Q'{t — T^ax)'^ and the X-ray 
data are consistent with r,^ax(BJD) = 2447650 ± 1200 and 
h' = -8.5(4.4) X 10"* d"^ or P33 = 2.0(1.0) x lO"^'' d d"^ 
(In both cases, the quoted errors account for the error in the 
de Jager spin ephemeris.) Of the two options, the quadratic 
fit is preferred, since it provides a slightly better fit to the 
data and produces phase residuals that are smaller (^ 0.07 
cycles for the quadratic term compared to ^ 0.13 for the 
cubic term) during the epoch studied by de Jager (1978.5- 
1992.6). In either case, the frequency derivatives are nega- 
tive and the period derivatives are positive, so during the 
epoch covered by the ASCA, XMM, and Chandra X-ray ob- 
servations (1995.8-2005.7), the white dwarf in AE Aqr was 
braking (slightly) harder than described by the de Jager spin 
ephemeris. 
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Table 2. Best-fitting cosine parameters A + Bcos(0a 



4>o)- 



Satellite 



Xu 



A 
(counts s"-*-) 



B 

(counts s"-*-) 



B/A 



<t>0 
(rel. cycles) 



Date of pulse maximum 
[BJD(TT) - 2400000] 



ASCA 1.43 

XMM 1.66 

Chandra 1.09 



0.799 ± 0.005 
8.232 ±0.026 
0.313 ±0.002 



0.103 ±0.007 
0.783 ±0.037 
0.047 ± 0.003 



0.129 ± 0.009 
0.095 ± 0.004 
0.151 ±0.009 



0.037 ±0.011 
0.146 ±0.007 
0.237 ±0.010 



50004.7037529 ± 0.0000044 
52221.5720142 ± 0.0000029 
53613.2767058 ± 0.0000039 
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Figure 3. Observed minus calculated spin phases of AE Aqr 
as a function of time {filled circles with error bars), the la error 
envelope of the de Jager spin cphemeris (dotted curves), the cubic 
fit to the de Jager optical datum and the ASCA, XMM, and 
Chandra X-ray data {solid curve), and the quadratic fit to the 
X-ray data [dashed curve). The epoch studied by de Jager is 
indicated by the thick horizontal line. 
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Figure 2. (a) ASCA, (b) XMM, and (c) Chandra spin phase 
folded count rate light curves of AE Aqr (filled circles with error 
bars), best-fitting cosines (solid curves), and mean count rates A 
(dotted lines). Each panel is scaled to ±25 per cent of the mean 
count rate. 



4 DISCUSSION 

Taking advantage of the very precise de Jager optical white 
dwarf orbit and spin ephemerides; ASCA, XMM, and Chan- 
dra X-ray observations spread over 10 yrs; and a cumulative 
27 yr baseline, we have found that in recent years the white 
dwarf in AE Aqr is spinning down at a rate that is slightly 
faster than predicted by the de Jager spin ephemeris. At 



the present time, the observed period evolution is consis- 
tent with either a cubic term in the spin ephemeris with 
P33 = 3.46(56) X 10"^^ d"^ centered on T^ax or an addi- 
tional quadratic term with P33 = 2.0(1.0) x 10"^^ d d"^ 
centered on T4ax(BJD) = 2447650 ± 1200. We consider the 
implications of each option in turn. 

If the observed O — C residuals are due to a cubic 
term in the spin ephemeris, it is possible to make an in- 
teresting comparison between the AE Aqr white dwarf X- 
ray pulsar and neutron star radio pulsars. In a vacuum, a 
rotating star with a misaligned magnetic dipole loses ro- 
tational kinetic energy via magnetic-dipole radiation at a 
rate L^d = -intl = 2/i^ sin^ Sn^/Sc^, where /i = BR^ 
is the magnetic moment, B is the surface magnetic field 
strength, R is the stellar radius, and 9 is the angle be- 
tween the rotation and magnetic axes. It is convenient to 
express this as f2 = —KQ" with K = 2/i^ sin'^ 9/3Ic^ and 
the so-called braking index n — 3. Ignoring the fact that 
the white dwarf in AE Aqr is not isolated [the hght cylinder 
radius (ric = c/fi = 1.6 x 10^^ cm) is comparable to the 
binary separation (a ~ 2 x 10^^ cm), so the white dwarf 
magnetic field will drag against the secondary, the mag- 
netic field of the secondary, and the mass lost by the sec- 
ondary] , if one assumes that the observed spin-down is due 
solely to magnetic-dipole radiation losses, one obt ains that 
^ Ri 1 X 10^* G cm^ or B ~ 50 MG Jlkhsanovll998h . compa- 
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rable to that observed in polars. However, one also obtains 
fl = ntf/n = 1.07 X lO"^'^ d"^ or P = -8.40 x 10"^" d'^, 
which resuhs is a decrease in O — C by 0.1 cycles in 10^ 
yrs, while we have observed an increase in O — C by 0.24 
cycles in 27 yrs. Expressed another way, we have derived not 
n = 3 but n ~ —42000. Clearly, the enhanced spin-down of 
the white dwarf in AE Aqr has nothing to do with magnetic- 
dipole radiation losses. On the other hand, if the observed 
O — C residuals are due to an additional quadratic term in 
the spin ephemeris, the spin-down rate is a modest 3.5 ± 1.8 
per cent greater than the rate derived by de Jager, which 
easily could be accommodated by, e.g., a small increase in 
the mass transfer rate from the secondary, leading to an 
enhanced spin-down torque on the white dwarf. 

While we have found that additional low-order terms 
in the Taylor series expansion of the white dwarf spin 
ephemeris adequately fit the trend of the O — C residuals 
of the ASCA, XMM, and Chandra X-ray light curves of AE 
Aqr, providing independent evidence of the nominal validity 
of the de Jager ephemerides, we caution that our results are 
based on only three data points spread over 10 yrs. First, it 
is possible that the spin evolution of AE Aqr is more com- 
plex than assumed, that O — C varies on shorter timescales 
and manifests larger excursions than sampled, and that the 
good quadratic and cubic fits to the O — C residuals of the 
X-ray data are fortuitous. Second, the quadratic and cubic 
spin ephemerides of §3 should be used with caution, since 
they have yet to be shown to have any predictive capability. 
Regular monitoring, in the optical, ultraviolet, and/or X- 
rays, is required to track the evolution of the AE Aqr pulse 
period, to determine if the white dwarf continues to spin 
down at the current rate, or if (as seems likely) it varies 
stochastically in response to changes in the mass-transfer 
rate, the varying efficiency of the magnetic propeller (e.g., 
the fraction of the mass-transfer rate that is expelled from 
the binary versus that accreted by the white dwarf), etc. 
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